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PACS 82 . 35 . Pq - Biopolymers. 

PACS 87. 16. Ka - Filaments, microtubules, their networks, and supramolecular assemblies. 

Abstract - Helically coiled filaments are a frequent motif in nature. In situations commonly 
encountered in experiments coiled helices are squeezed flat onto two dimensional surfaces. Under 
such 2-D confinement helices form "squeelices" - peculiar squeezed conformations often resembling 
looped waves, spirals or circles. Using theory and Monte-Carlo simulations we illuminate here the 
mechanics and the unusual statistical mechanics of confined helices and show that their fluctu- 
ations can be understood in terms of moving and interacting discrete particle-like entities - the 
"twist-kinks". We show that confined filaments can thermally switch between discrete topological 
twist quantized states, with some of the states exhibiting dramatically enhanced circularization 
probability while others displaying surprising hyperflexibility. 



Introduction. — Helically coiled filaments are found 
everywhere in living nature. The list of examples is close 
to innumerable with the most prominent ones: FtsZ pQ, 
Mrb [2], bacterial flagella [3j|4], tropomyosin [5] and in- 
termediate filaments [6]. More recently microtubules 
were suggested to spontaneously form large scale super- 
helices [7J|S]. Even whole microorganisms exhibit helic- 
ity inherited from their constituent filaments [S]. The 
superhelicity of filaments is in some cases of strong evo- 
lutionary benefit as in the example of swimming bacte- 
ria utilizing the rotational motion of their helical flagel- 
lar filament for propulsion [ 1 OJ and tropomyosin's helical 
" Gestalt-binding" around actin [5]. In other cases, like 
for microtubules the purpose of superhelicity remains so 
far unknown [8] . Paralleling biological evolution artificial, 
man made helically coiled structures have been created 
including coiled carbon nanotubes [IT], DNA nanotubes 
[T2"] . coiled helical organic micelles 

As happens often in experiments our justified desire to 
simplify observation conditions by confining filaments to 
a surface (coinciding with the focal plane) changes the 
physical properties of the underlying objects in initially 
unanticipated but physically rather interesting manner. 
With filament helices being such a profoundly ubiquitous 
structure it is the purpose of this work to investigate the 
rich physical effects of helical filaments confinement. As 
we will show, the confinement changes dramatically the 
shape as well as statistical mechanics of the confined he- 



lix generating several notable and surprising effects: a) 
Enhancement of cyclisation probability, b) Enhancement 
of end-to-end fluctuations and c) Generation of confor- 
mational multistability (despite apparent linearity of con- 
stitutive relations). We will see that the conformational 
dynamics of confined helices is most naturally described in 
terms of discrete particle like entities - the " twist kinks" , 
cf. Fig. 1. We show that these "twist kinks" are com- 
pletely analogous to overdamped Sine-Gordon-kinks from 
soliton physics [14] as well as loops in stretched elastic 
filaments [15] . These analogies will help us to rather intu- 
itively develop a phenomenological understanding of the 
underlying physics. 

The Phenomenology of Squeezed Helices. 

Confined biofilaments throughout literature exhibit often 
abnormal, wavy, spiral, and circular shapes that appear 
not be rationalized by the conventional Worm Like Chain 
model. This riddle of peculiar filament shapes is the start- 
ing point of our investigation. In this letter, we propose 
a new augmented model of confined intrinsically curved 
and twisted chains that leads to a variety of 2-D shapes 
matching experimental observations. The corresponding 
2-D geometrical curves we call squeezed helices or more 
briefly - squeelices. Filaments are modelled as 2-D con- 
fined Helical Worm Like Chains (cHWLC) with bending 
modulus B, twist modulus C and a preferred curvature o>i 
and twist L03. 
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Phcnomcnologically (as detailed further below) the 
main effect of confinement is to introduce narrow regions 
where the twist is highly concentrated and the curvature 
flips, cf. Fig. 1. These curvature flip points we will call 
" twist-kinks" in reference to the concept of kinks in soliton 
physics [14j . Depending on the control parameter 7 ~ 
two regimes can be distinguished: 

i) 7 > 1 : Twist kinks having a positive self-energy are 
essentially expelled and can only be thermally activated. 

ii) 7 < 1 : Twist kinks have a negative self-energy and 
the ground state involves a finite density of twist kinks. 

In both cases the generic shape motive is a circular arc 
section or in special cases closed circles. If present, twist 
kinks separate arcs of opposite curvature orientations re- 
sulting in a wavy undulatory 2-D shape. For a squeelix 
with 7 > I, twist expulsion results in a circular arc shaped 
ground state -a feature that can favor the occurrence of 
closed filaments, which are very unlikely in the free helical 
state (in absence of confinement). 

Our study is based on the analytical analysis of the 
confined helical WLC Hamiltonian and numerical Monte 
Carlo simulations using the density of states method [IS] . 
The simulations, primarily designed to illustrate circular- 
ization enhancement, focus on the most curious case of 
twist expulsion where the ground state becomes circular 
up to (weak and localized) edge effects. In this regime, 
besides the ground state, the chain will also comprise ex- 
cited (wavy looking) states involving a discrete number of 
thermally activated twist-kinks. Because an isolated twist- 
kink almost freely diffuses along the filament an excited 
state can be " hyperflexible" meaning that large shape fluc- 
tuations are induced by simple translational sliding motion 
of the twist-kink along the contour. For sets of parameters 
not resulting in twist expulsion (7 < I), the ground state 
itself is periodically wavy. This case where a squeezed he- 
lix forms a wavy periodic 2-D structure appears less excit- 
ing as effects like circularization and hypcr-flcxibility will 
be absent. For this reason we will entirely focus in this 
very first study on the more illuminating and surprising 
limit 7 > 1. 

A note for the experimentalist concerning non-ideal fila- 
ments (in real world experiments) seems appropriate here. 
Bacterial flagella and microtubules are known to be helical 
filaments, but their mechanical characteristics are gener- 
ally (in experimental practice) not quite uniform along the 
filament. There can be frozen-in lattice defects like a jump 
in the number of protofilaments for microtubules [IT] or 
annealed defects like a boundary between different poly- 
morphic states in bacterial flagella. In the cases of interest 
the mechanical parameters are uniform within some cor- 
relation length, which may even be larger for cleanly pre- 
pared filaments than the contour length under scrutiny. 
Our model is generic in the sense that it applies to a uni- 
form sub-filament. A clean observation and comparison 
of the shapes described below further supposes that shape 
equilibration is possible with uniform mechanical param- 



eters. 
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Fig. 1: Squeezed helical filaments (" squeelices" ) are frequently 
observed under experimental conditions, a) Polymorphic he- 
lical bacterial flagella emerging from a bacterium [18] (Scale 
bar 1.3/j.m). b) A coiled carbon nanotube when adsorbed onto 
the substrate from [TT], (Scale bar 0.5/xm) c) Helical DNA nan- 
otubes [12] (Scale bar, 0.5/^m). c) In simulations squeelices ap- 
pear in coexisting wavy-circular shapes and display discrete 
"twist kinks" (curvature inversion points) indicated by arrows. 
Twist kinks - discrete non-linear elementary excitations of con- 
fined helices - are characterized by a rapid inversion of the nor- 
mal vector along the contour and by a localized twist-angle ip 
variation 



The Squeezed Helical Chain Model: Free vs 
Confined. — The variety of squeelical shapes in Figs. 
1 and [5] can be understood by simple considerations about 
their energy. The clastic energy of a helical wormlike chain 
can be written as a function of the local curvatures fii^ 
and twist O3 as: 



2 



L/2 



b [(fi! ~ u^) 2 + n 2 2 ] +c(n 3 -LJ3) 2 ds, (1) 



-L/2 

The 3-D ground state is a helix of radius R 
and pitch H = satisfying the preferred curvature 

and twist everywhere. To proceed and access the fila- 
ment shapes, it is convenient to express the VLi through 
the Eulcr angles Q,\ = <p' sin 9 sin ip + 9' cosip, SI2 = 
<p>' sin 9 cos ip — 9' sin-0 and 17 3 = <p' cos 9 + ip'. We con- 
strain the chain to a plane by imposing 9 = tt/2. The 
elastic energy of such a confined squeelix can be recast as: 



E = 



D 



L/2 



(V) 2 - 2wi^ sin i> + lo 2 + c (if/ - 0J 3 f) 



da 



-L/2 



(2) 

where c = C / B and s is the arc length parameter along the 
chain. Now cf)' 2 = Q\ + fl 2 , = k 2 gives the local curvature 
of the chain in the plane and ip' = is its local twist. 
Interestingly it is seen from Eq[5] that the curvature and 
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the twist are now coupled through oj\ - which turns out to 
lie at the heart of most phenomena discussed here. The 
optimal shape satisfies the Euler-Lagrange equations: 



4> = u>i simp. 

ib" + ^sia2ip = 0. 
2c 



(3) 
(4) 



If no external torque is applied on the chain, the twist must 
satisfy the boundary conditions ip'(—L/2) = ip'(L/2) = 
L03. A quick glance at Eq. |4] reveals that ip(s) is the solu- 
tion of a pendulum equation and thus is a Jacobi elliptic 
function [T5] . Eq. [3] shows that under confinement the 
curvature becomes "slaved" by the twist and is no more 
an independent parameter as for a free chain. This obvi- 
ously is at the origin of the localization of twist. For the 
sake of simplicity, instead of solving Eqs. 13141 directly, we 
will gain more physical insight by rewriting Eq. [5] in terms 
of a WLC under tension (WLC-T) for which it is easy to 
develop intuition. For this we plug in Eqj3]into the energy 
Eq[2]and introduce the new angle •& = 2ip — it to get: 

L/2 

E{$)= J Qi?9' 2 + F(l-costf)^ ds 

12 

[0 (L/2) -0 {-L/2)], (5) 



-L/2 

M 
~2 



up to an inessential constant term. Here we introduced 
A = C/4 as the effective bending modulus of the WLC-T 
and F = Buj\/A as the effective external force acting on 
it. The integral term is precisely the WLC-T energy. The 
last term with M = u^C represents the torque exerted on 
both ends of the chain. If large enough, this term enforces 
n extra trapped turns/loops •& {L/2) — d (—L/2) = 2nn or 
in the ip representation, n twist-kinks along the squeelix. 

The phenomenology of the stretched chain is com- 
mon knowledge |15j . which makes the mapping attractive. 
The decay length of any localized distortion (correlation 
length) A is the loop size of a loop grown against the ten- 
sion F (in the loop picture) or kinks-spatial extension (in 
the twist-kink picture): 



A 




(6) 



The WLC-T loop stores a typical energy ~ XF. If the 
work of the external torques 2ttM per loop reduces suffi- 
ciently the loop energy, loops form spontaneously (other- 
wise they can be thermally activated). We hence conclude 
that the ground state is wavy for M > XF, which trans- 
lates into I > r in the cHWLC language, and circular 
otherwise. 

An overestimate of the energy stored in the loop 
2\/2tx\/ AF is obtained assuming a circular loop. The 
WLC-T featuring loops in the long chain limit has been 
studied in detail in rcf. 15 from which we borrow the more 



precise expression Ei oop = 8v AF. After subtracting the 
work of the torques we arrive at the effective self-energy 
of a single twist-kink: 



Eikink = Wl~ 1)ttCoj 3 



where 



7 



(7) 



(8) 



is the twist-kink expulsion parameter introduced previ- 
ously. The shape of an isolated twist-kink is easily ob- 
tained upon integrating Eq[?l for a twist kink localized 
around s = 0: 

s s 
simp (s) = tanh — and <p> (s) = U)\ tanh — (9) 
A A 

a result we could also have transposed from [15]. From 
Eq|9] it is clear that provided L > J the twist kink is 
localized and separates two curvature flipped regions of 
almost constant curvature oj\ |20j . 

Let us now turn to the main focus of this paper, the case 
of twist expulsion (7 > 1), which maps to the ground state 
to a WLC-T without loops. In a consistent thermody- 
namic picture the twist-kink density should be calculated 
including shape fluctuations, in particular the coupling of 
the twist kink to small thermal deformations. These linear 
fluctuations -which in the soliton picture are often referred 
to as "phonons" [14] - are expected to only moderately 
renormalize the kinks free energy (by ~ ksT per kink). 
Fortunately intuitive reasoning and simple geometry al- 
lows us to calculate filament's shape fluctuations associ- 
ated with the free sliding motion of the narrow twist-kink 
along the filament in the excited state. The associated 
soft mode (kink's position) is at the origin of the surpris- 
ingly large fluctuations in the excited state reported in 
simulations below. 

The opposite regime 7 < 1 where the squeelix shapes 
become wavy can be described along similar lines with 
the twist kink density being limited by mutual kink- 
kink repulsion. Generally speaking for 7 < 1 the gas 
of twist kinks seizes to be ideal, giving rise to reduced 
compressibility and in turn weaker extension fluctuations 
of the squeelix. It can be shown that the pair repul- 
sion between twist-kinks decreases with their distance d 
as Uint ~ 7rCo;3 % /7/(d/A) with f(x) ~ l/x for i«l and 
f(x) ~ e~ x / 2 for i> 1. The dense twist kink regime de- 
serves special consideration, a more detailed description 
explicitly involving Jacobi functions will be given else- 
where. 

Simulation. — Equipped with this intuitive picture of 
a squeezed chain as a collection of weakly interacting ID 
" particles" (twist-kinks) whose positions along the contour 
relate to chain's deformations in simple manner (Eq. |9| 
we move on to investigate the formulated hypotheses with 
a Wang-Landau type Monte-Carlo simulation [16]. We 
model the system as a discrete helical WLC consisting of 
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iV-monomers [H] subjected to a discrete version of the 
Hamiltonian given in EqUJ The chain is further confined 
by a harmonic potential so that each monomer located 
at distance z away from the z = surface experiences 
potential E con f = Kz 2 with K = 25k B T/b 2 , with the 
bond length set as b = 1. For K = the confinement 
vanishes and we recover the free chain statistics in 3-D. 

The main output of the simulation is the density of 
states (DOS) [16]. The equilibrium conformational statis- 
tics can be obtained from the DOS after Boltzmann 
weighting according to the Hamiltonian. In addition to 
the energy E, the joint DOS g(E, D, n) was sampled with 
respect to the end-to-end distance D and to the "number 
of kinks" n, where n = ~ 2 fia^. The latter is indeed 
a statistical analogue of the discrete twist kink number 
in the theoretical consideration above. Sampling of two- 
dimensional or higher histograms is often computationally 
more demanding than one-dimensional histograms. For ef- 
ficient sampling, we used the global update method intro- 
duced by [27] . The simulations were performed for several 
chain lengths, N = 16, 32, 64 and with fixed mechanical 
parameters B = 50, C = 25 6fcsT throughout. The essen- 
tial information on the confined cHWLC gained from our 
simulations are summarized in free energy maps F(n, D) 
(see figgpO). 

Results. — Figure [2] demonstrates the typical changes 
of chain's conformation triggered by confinement. For val- 
ues 7 > 1, the confined helix assumes an approximately 
circular shape as shown in Fig[2](a) . The binormal vectors 
pointing all in the same direction (in red) indicate that 
twist is expelled from the chain. The normal vectors are 
in-plane and pointing to the center of circle (in blue). For 
7 < 1, on the other hand, the chain assumes a stretched 
and twisted undulatory wavy shape as seen in Fig[2](b). 
The wavy shape consists of a discrete number of localized 
twist kinks as anticipated from theory. At the location of 
twist-kinks the in-plane normal vectors and the in plane 
curvature </>' flip their signs. The curvature between the 
twist kinks is approximately constant w uj\ as expected 
from the theoretical shape of single twist kinks Eq. [9] 

Hypercyclisation. The phenomenological theoretical 
arguments above predict that a significant circularization 
enhancement should be observable for a twist-expelling 
chain (7 > 1). We subject this hypothesis to a simulation 
test. For the geometric parameters ojiN = 4.8 and 
0J3N — 3.2 and elastic constant ratio B/C = 2 the 
twist expulsion parameter is well over unity 7 = 1.82. 
In absence of confinement, a chain with these param- 
eters assumes a helical shape with a radius given by 
R = 0.14A6 and a pitch H = 0.62 Nb. When on the other 
hand such a chain becomes strongly confined to 2D , it 
morphs into a slightly open but approximately circular 
shape with an end-to-end distance approximately given 

b y m ~ z^Tv s * n (^) = °- 28 - In such a case ' the 

newly formed ground state under confinement seems to 




Fig. 2: Formation of two different classes of conformations of a 
confined helix depending on the twist expulsion parameter 7. 
The contour is represented by a tube with a normal (in red) 
and a binormal vector (in blue), (a) and (b) typical helical 
chains in 3-D before confinement. When confined, (a ) un- 
twisted circularized shapes (7 > 1) and (b ) twisted and wavy 
shapes (7 < 1) emerge. 



promote circularization. To quantify the circularization 
enhancement we define the circularization probability 
PP(D cap ) given by the likelihood to find both chain ends 
within a fixed (small) capture distance D cap 



cap } 



Pd(D)dD. 



(10) 



where p d (D) = ±ff g (E, D 7 n) e~ E / kBT dEdn , Z = 
J J fg(E,D, n) e~ E t kBT dEdndD is the radial probabil- 
ity density for a given end-distance D and the index d = 2 
or 3 stands for the dimensionality of the chains embed- 
ding (confined in 2D or free in 3D). FigO shows the radial 
distribution functions p 2 (D) (2-D) and p 3 (D) (3-D) at var- 
ious temperatures. The maximum of p%(D) (blue line) is 
located at about D = 0.55Nb close to the expected helical 
pitch. In contrast the largest population of the confined 
chains is located around the peak of P2(D) (blue line) at 
Do/Nb « 0.28 which corresponds to the expected 2-d cir- 
cular ground state. At higher temperature, chain becomes 
more flexible and the end-to-end distribution broadens due 
to the thermal fluctuations. 

Integrating the probability density EqllOl we can es- 
timate the circularisation enhancement factor /O = 



P% (D C ap)/ P-$ (D ca p) when the chain is squeezed to 2D 
as a function of the chains bending persistence length 
lp = B/ksT. As suggested on theoretical grounds the 
closing probability is greatly enhanced due to the confine- 
ment at lower temperatures (or increased chain stiffness) 
- with fO well in excess of 1000 (Figdfc). 

Multistability and Hyperflexibilty. The theoretical 
considerations of multi-kink solutions presented above in- 
dicate that confining a helical chain should generate a 
complex energy landscape with many coexisting discrete 
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Fig. 3: Confinement enhances circularization. The end-to- 
end distance D distribution as a function of chain's bending 
persistence length l p — B/kBT (a) in 2D (a) and (b) 3D 
(N = 32b, wi = 0.156 -1 , oj-a = 0.1b" 1 , and 7 = 1.82 > 1). The 
enhancement of circularization probability upon confinement 



/ u = Pr^ (D cap ) /Pgj (D cap ) as a function of the normalized 
chain length Nb/l p with capture distance D cap — 0.5b. 



states. Some of these states, in particular those comprising 
only a few kinks, should exhibit anomalous hyperflexible 
behavior due to the energetically cheap displacement of 
kinks. For instance the motion of a single kink on a chain 
of length L and curvature ui\ satisfying a close to "reso- 
nance" condition oj\ w fc^S (k = 1,2... ) should give rise 
to most dramatic end fluctuation effects. 

To test this we calculate the joint DOS of three vari- 
ables g(E,D,n) and the free energy landscape given by 
PF[D,n) = -ln[J2 E g(E,D,n)e-^ E ] (Fig. 4a,b). For 
parameters u>\ = 0.266~ 1 and CJ3 = 0.16 -1 , in 3D (con- 
finement free case) we would obtain a helix with helical 
radius R = 8.096 and pitch H — 3.356. In contrast to 
that as expected from the large twist expulsion param- 
eter 7 = 5.47 the chain's 2D ground state shape under 
confinement becomes a coiled double circle with roughly 
two turns. Besides the ground state, the free-energy as a 
function of n indicates the existence of further ("excited") 
low energy states with n = 1, 2, ... . (Fig. [4j). These states 
with distinct n are separated by small free energy barriers. 
The first excited state n = 1 for instance has a free energy 
difference of only AGikink ~ 4fcsT a figure that is not 
far from the theoretical estimate of the kink free energy 
AGifejnfe = E lkink - k B T\nN « 6M B T where the In JV 
(N=48) term accounts for kinks positional entropy gain 
along the discrete positions of the chain [52] . 

The n = 1 state (cf. Fig. 2{c)) exhibits an enhanced 
end-end distance fluctuation ((SD) 2 ) in phenomenologi- 
cal agreement with the mobile kink interpretation. For 
the given length (L = 486) and curvature we expect 
an almost flat free energy landscape as a function func- 
tion the end-end distance D in the range D £ [0,D max ] . 
Here D max ~ 4wJ~ 1 = 15. 46 is the maximal extension for 
the n — 1 state which is reached when the twist-kink is lo- 
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Fig. 4: (a) The free energy landscape of a fluctuating squeelix 
as a function of the end-to-end distance D and the number 
of twist kinks n (N = 486 , wi = 0.266 -1 , oj 3 = 0.16" 1 , and 
7 = 5.47 > 1). The ground state E g (n = 0) is a two-turn cir- 
cular shape without kinks. The excited states n > 1 comprise 
one or more very mobile kinks giving rose to strong shape fluc- 
tuations with states E\ a and Eib having identical energies, (d) 
Contour plot of free energy for 2-d chain with indicated values 
free energy values in units of ksT. (c) Mean of end-to-end 
distance (D) and its standard deviation {(SD) 2 ) as a function 
of n (dual axis graph). 



cated at a position w L/4 from any of the borders- as seen 
from straight forward geometric reasoning. The standard 
deviation expected in this case ((SD) 2 ) n=1 rj 286 2 com- 
pares well with the simulation result ((SD) 2 ) n= i rj 266 2 . 

The higher states n > 2 display again lower fluctua- 
tions. This is in agreement with the interpretation that 
with growing kink density their repulsion and eventually 
mutual confinement become important. 

Conclusion. — We have shown that the conceptu- 
ally simple procedure of planar confinement transforms 
a simple mundane object - a helically coiled filament- 
into a complex multistablc, anomalously fluctuating fila- 
ment. The statistical mechanics of this exotic object -the 
"squeelix"- can be qualitatively understood in terms of 
the motion of discrete particle-like entities corresponding 
to sharp curvature inversion points - called " twist kinks" . 
At low twist kink concentrations they move almost freely 
along the chain and induce anomalously strong conforma- 
tional fluctuations- notably deviating from wormlike chain 
behavior. The "squeelical" shapes formed under confine- 
ment range from almost ideally circular to wavy depending 
on the value of a single dimensionless " twist expulsion" pa- 
rameter 7 = ^X 1 ,^ ■ The latter parameter depends both on 



the chain's elastic moduli (flcxural modulus B, twist mod- 
ulus C) and the geometric properties (intrinsic curvature 
u>\ and twist W3). For 7 > 1, the twist is curiously expelled 
from the chain. Under these conditions the squeelix be- 
comes almost circular (up to minor end-effects) and the 
circularization probability can be dramatically enhanced. 
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In the other limit 7 < 1 the squeelix comprises densely 
packed twist-kinks in its ground state which is now wavy. 

We may speculate whether abnormal wavy shapes or 
enhanced closure of confined filaments found in literature 
are fingerprints of an otherwise hidden helical superstruc- 
tures whose microscopic origin should be elucidated for 
each type of filament, case by case. Such peculiar behav- 
ior under confinement is observed for essential filaments 
like microtubules j7], F-actin [28] and possibly for inter- 
mediate filaments [31J . It has not escaped our attention 
that actin filaments circularize on stunningly small scales 
(~5/Ltm length rings) [28], |29j and exhibit wavy periodic 
tangent-correlation functions in narrow, flat -channels (cf. 
Fig .6 in |30| ) - both phenomena that neither can be un- 
derstood within the naive WLC model. Similarly inter- 
mediate filaments under confinement show wiggly periodic 
shapes [31] suspiciously reminiscent of squeezed helices. 

It is our feeling that the here illustrated phenomena of 
multistability, hyperflexibility and enhanced chain circu- 
larization are just the tip of an ice-berg and that floppy 
"squeelix" concept will help us to unravel a number of 
mysteries from biofilament realm hidden out in literature. 
We have evidence that additional lateral confinement of 
filaments in narrow channels, as encountered in microflu- 
idic devices |3I] , further enhances the visibly wavy shapes 
for helical filaments and gives easier access to the under- 
lying mechanical parameters. The potential of the 2-D 
and especially double confinement experiments (in mi- 
crochannels) has been vastly underestimated so far. The 
"squeelix" phenomenology laid out here will serve us as 
a "dictionary" to decode these peculiar observations in 
forthcoming works. 
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